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1.      INTRODUCTION 

This  report  is  a  presentation  of  one  part  of  a  study  on  the  contaminating 
backflow  from  the  exhaust  plume  of  a  large  space-based  HF  laser  (Figure  1). 
The  exhaust  plume  is  an  underexpanded  supersonic  ring-jet,  designed  to  stay 
clear  of  the  spacecraft  by  maintaining  a  Prandtl-Meyer  turning  angle  at  the 
nozzle  lips  of  well  below  90°.  However,  it  is  well  known  from  experience  with 
rocket  plumes  in  space [1.2]  that  cavity  regions  (where  continuum  gasdynamic 
theory  predicts  vacuum)  are  filled  with  a  free-molecular  flow.  This  back  flow 
is  largely  due  to  viscous  effects,  which  give  rise  to  a  "spill-over"  of  the 
boundary  layer  around  the  nozzle  lipt^].  Assuming  the  boundary  layer  can  be 
eliminated  (e.g.,  an  expanding  step  design  of  the  nozzle  lip),  there  are 
two  more  mechanisms  which  lead  to  backflow:  thermal  backscattering  and 
scattering  by  ambient  molecules  traveling  at  orbital  speeds.  Since  these 
effects  are  a  small  perturbation  to  the  exhaust  flow  field,  they  can  be 
considered  independently  (the  total  backflow  will  be  a  superposition  of 
contributions  due  to  these  two  effects).  As  a  first  phase  of  our  broader 
study,  we  consider  solely  the  contribution  of  thermal  backscattering  to  back- 
flow  from  a  ring-plume  of  an  HF  laser,  via  a  simple  model  of  molecular 
effusion  from  a  breakdown  surface,  fashioned  after  ideas  suggested  by 
Nollert^i.  Our  results  indicate  that  the  backflow  of  the  heavier  contaminants 
(HF,    DF)    due   to   thermal   backscattering   is    negligible. 

Naturally,  our  study  pertains  to  presumably  typical  operating  conditions 
of  the  HF  laser.  These  operating  conditions  were  largely  determined  from  a 
report  on  some  HF  laser  tests  conducted  at  TRW  in  1971 [5]  (in  particular, 
Table    5,    Test    III,    of   this    report).      The   typical   parameters    at    the    nozzle'exit 


are 


[H]  =  .091,  [HF]  =  .091, 
[H2]  =  .104,  [DF]  =  .135, 
[He]  =  .579 

T  =  1.54  (assuming  ideal  gas) 

Mi  =  4.0 

WA  =  7.27  [kg/kg  mole]   .  (1.1) 

T0  =  2300  [K] 

P0  =  0.0075  [kg/m3] 

D  =  2.5  x  10"10  [m] 


Composition  (mole  fractions): 

Specific  Heats  Ratio: 

Mach  Number 

Average  Molecular  Weight 

Stagnation  Temperature 

Stagnation  Density 

Molecular  Diameter,  assuming  it  is 
uniform  for  all  species  (hard- 
sphere  collisions) 

The  exit  Mach  number  can  be  chosen  higher  than  Mi  ■  4,  but  not 
considerably  lower  than  this  value,  since  Mj  *  4  results  in  a  modest  clearance 
angle  of  41°  between  the  limiting  (vacuum)  characteristic  of  the  lip-centered 
rarefaction  fan  and  the  spacecraft.  We  assume  isentropic  flow  throughout  the 
diffuser^J,  so  that  upon  specifying  the  composition  and  flow  variables  at  the 
diffuser  inlet,  along  with  Mi  at  the  diffuser  exit,  the  exit  flow  is  fully 
determined.  One  exception  to  this  definition,  however,  is  the  stagnation 
temperature,  which  was  estimated  as  T0  =  1400  [K]  at  the  diffuser  inlet i^]. 
We  set  T0  =  2300  [K],  which  corresponds  to  complete  hydrogen  recombination, 
even  though  the  flow  in  the  diffuser  is  of  a  nearly  frozen  composition  due  to 
the  low  rate  of  hydrogen  recombination  1 5] .  The  reason  for  this  choice  is  that 
given  the  uncertainty  in  determining  T0,  which  results  from  an  uncertainty  in 
the  degree  of  hydrogen  recombination,  it  is  the  most  pessimistic  choice, 
resulting  in  higher  thermally  backscattered  flux. 


The  model  that  we   propose  for  evaluating  the  backscattered   flux  arriving 


at  the  spacecraft  (Figures  1,  2)  is  based  on  the  effusive  breakdown  surface 
concept  suggested  by  Noller^J.  The  gradual  transition  from  continuum  to 
collisionless  flow,  which  invariably  takes  place  at  the  outer  fringes  of 
exhaust  plumes  having  a  near-vacuum  background  environment,  is  replaced  by  an 
abrupt  change.  We  assume  that  the  flow  regime  in  each  stream  tube  changes 
from  continuum  (with  local  thermodynnamic  equilibrium)  to  collisionless,  upon 
crossing  some  breakdown  surface. 

An  important  simplification  is  introduced  in  the  case  of  a  large-radius 
spacecraft  (about  2.5[m]),  by  observing  that  the  temperature  along  the 
breakdown  surface  decreases  so-  sharply  with  the  distance  from  the  .  nozzle  lip, 
that  the  segment  contributing  significantly  to  thermal  backscattering  is  only 
about  0.01  to  0.1 [ra]  long.  Consequently,  the  lip-centered  rarefaction 
ring-fan  may  well  be  approximated  by  the  standard  (planar)  Prandtl-Meyer 
flow  field. 

The  structure  of  this  report  is  as  follows.  The  breakdown  surface  and  the 
molecular  effusion  flux  from  it  are  obtained  in  closed-form  expressions  in 
Section  2.  Section  3  is  devoted  to  the  spatial  integration  scheme,  which  is 
the  evaluation  of  the  flux  arriving  at  a  certain  point  on  the  spacecraft. 
Results  of  flux  distribution  along  the  spacecraft  for  the  presumed  laser 
operating  range  are  presented  and  discussed  in  Section  A,  followed  by  a 
critical  examination  of  the  breakdown  surface  model.  Conclusions  are  given 
in  Section  5,  and  Section  6  is  a  list  of  references.  The  code  RINGBD,  which 
computes  the  flux  by  numerical  integration  over  the  breakdown  (effusing) 
surface,  is  given  in  Appendix  A. 


2.   BREAKDOWN  SURFACE  AND  EFFUSION  FLUX 

Our  model  for  the  thermally  backscattered  flux  arriving  at  the  surface  of 
the  spacecraft  is  essentially  a  modification  of  Noller's  concept  of  a 
breakdown  effusive  surface^.  We  substitute  his  definition  of  a  breakdown 
surface  by  a  similar  one  introduced  by  Bird^6'  Section  8.3] #  we  obtain  the 
one-sided  effusion  flux  from  the  breakdown  surface  by  integrating  over 
velocity  space  as  suggested  by  Nollerl^J,  except  for  the  fact  that  we  compute 
flux  rather  than  density  and  we  also  consider  the  flux  of  species  having 
molecular  weight  different  from  the  average.  In  the  following,  each  one  of 
these  steps  is  described  in  some  detail,  beginning  with  the  breakdown  surface. 

As  mentioned  in  the  introduction,  the  lip-centered  rarefaction  fan  is 
approximated  by  a  planar  Prandtl-Meyer  flow  field  (Figure  2).  The  standard 
expressions  for  this  flow  field  have  Mach  number  (M)  as  the  independent 
parameter,  thus  M  varies  between  M  =  M^  at  the  exit  and  B  ♦  •  at  the  limiting 
(vacuum)  characteristic.  (Index  1  always  refers  to  exit  conditions,  i.e.,  to 
parameters  evaluated  at  M  -   M  ). 

ij/(M)  =  -  c(M)  +  5l  +  v+   J 

C(M)    =    rW2   ARCTAN[r-1/2(M2-l)1/2]  H^) 

u(M)   -  ARCSIN(M-1)  (2,1) 

8(M)  =  ^(M)  -  u(M)  ' 
where  \\>     is   the  angle   of  characteristic   lines,  and  9   is  the  angle   of  the 
velocity  vector  (or  streamline). 


Adopting  Bird's  definition  of  a  breakdown  parameter,  which  was  first 
introduced  in  conjunction  with  a  spherical  source  flow^°  »  Sectlon  8.3]  and 
later  was  shown  to  be  meaningful  also  in  a  Prandtl-Meyer  flow^'J,  we  defi 
the   breakdown   surface      as    having   a   constant   value      of    B,    where    B   is      given   by 


ne 


=  U  j_  Idp 
v   p  I  dS 


(2.2) 


Here   p,  U,   v,  S  are  local  flow  density,  speed,   collision   frequency,  and 

coordinate  along  streamlines   (thus   restricting   this   definition  of   B  to 

stationary  flows).    From  the  geometrical  relationships  in  a  Prandtl-Meyer  fan 
(Figure  2)  and  from  (2.1)  we  get: 

§  -  -  (1/R)  {&)    sin  u  -  "  -^   [M-1(M2  -  1)1/2]  (P/R) 

1  (2.3) 

P(M)  =  Po  U  +^M2)   Y_1 

Using    the   expression   for    collision   frequency '"  -I  : 

v0  =   4(rr/Y)1/2      (N0D2C0)  (2.4) 

where   N0,C0,D   are   stagnation   number   density,    stagnation   sound   speed,   molecular 
diameter,    and   using   U   =  MC    in   conjunction   with    (2.2)    and    (2.3),    we   get: 

_1 


RB(M)    =   (BN0D2)"1  tyzLr     [M2-l)1/2    (1+JCi  M2)   y-1  (2.5) 


This  expression  is  almost  identical   to  that   of  Birdt?],  the  main  difference 
being   in  assuming   a  constant  collision  diameter  (hard   spheres),   which  we 


believe  to  be  commensurate  with  the  overall  crucieness  of  the  model.  The 
breakdown  surface  as  defined  by  equations  2.5  and  2.1,  starts  at  point 
[Rg(M  ),  i|»(H  )]  on  the  exit  characteristic  M  =  M  ,  which  we  refer  to  as  the 
intial  point  (see  Figure  4).  However,  a  breakdown  in  continuum  flow  also 
takes  place  on  the  segment  of  the  exit  characteristic  between  the  corner  and 
the  intial  point,  since  the  value  of  the  breakdown  parameter  there  (Equation 
2.2)  is  clearly  larger  than  the  value  of  B  used  in  defining  the  breakdown 
surface  (Equation  2.5).  Hence,  the  breakdown  surface  defined  by  (2.5)  should 
be  supplemented  by  that  segment.  We  refer  to  the  combined  surface  as  the 
augmented  breakdown  surface.  The  segment  on  the  exit  characteristic  is 
referred  to  as  the  supplementary  breakdown  surface. 

The  one-sided  directed  effusion  flux  is  defined  as  the  number  flux  of 
molecules  per  unit  area  of  an  area  element  normal  to  the  flux  direction,  per 
unit  solid  angle  about  the  flux  direction.  It  is  obtained  as  a  function  of 
local  llach  number  and  the  angle  <  between  the  flux  direction  and  the  local 
velocity  vector,  by  repeating  Noller's  velocity  integration  scheme^'™'  (6)  J, 
with  an  added  factor  of  molecular  speed  in  order  to  obtain  flux  (rather  than 
density  as  in  Noller's  work).  The  resulting  expression  for  species  i  is 
readily  obtained  by  using  standard  definite  integrals: 


Fi<M)  =  ^  NQC0  (WA/Wi)l/2  (1  +121  M2)   2(Y_1) 

[(2ytt3)"1/2  (l  +  (l/2)yM2cos2<)  EXP  (-  yM2/2)  + 

(2.6) 

(2tt)_1  M  (3/2  +  (1/2)yM2cos2<)  cosk  EXP  (-  (l/2)yM2  sin2<) 

1/2  ~ 
ERFC  (-  (Y/2)     M  coskJ] 

00 

ERFC(v)    =   2   tt_1/2    J  exP    (-x2)dx  (complementary   error   function) 


5  =  (W±/WA)1/2  M 
hj  -  Mole  fraction  of  species  i.  (2.6  Continued) 

W^  -  Molecular  weight  of  species  i. 

The  dependence  of  F^(M)  on  the  flux  angle  <  is  so  sensitive  that  for  some 
Mach  number  around  M  *  10,  the  backflow  (in  the  typical  operating  range)  is 
virtually  negligible.  In  the  following  section  we  describe  how  the  flux  F^(M) 
is  integrated  over  the  augmented  breakdown  surface,  yielding  the  backscattered 
flux  arriving  at  the  surface  of  the  spacecraft. 


3.   FLUX  INTEGRATION 

The  effusion  flux  Fj_(M)  given  by  (2.6)  above,  is  defined  in  such  a  way 
that  the  number  of  molecules  effusing  from  an  area  element  AAg  of  the  break- 
down surface  and  arriving  at  an  area  element  AAg  on  the  spacecraft  (per 
second),  is  given  by: 

Fj[(M)  (AAg  COS  aB)  (AAS  COSas)  L~2  [molecules  per  second]        (3.1) 

o  O 

where  ag,  ao  are  the  angles  between  the  line-of -sight  Lgc  (Figure  3)  and  the 
normals  to  the  breakdown  surface  and  the  spacecraft  surface  respectively.  Lgg 
is  the  distance  between  the  elements  AAg  and  AAg.  Dividing  equation  (3.1)  by 
AAg  and  integrating  over  the  breakdown  surface,  the  flux  per  unit  area  of 
the   spacecraft   is    given   by: 

Qi  =   /F1(M)COSaB  C0S<*s  L~2   dAB        [Molecules    per   second    per  m2]         (3.2) 

The      integration   scheme      for   Q^   over    the   breakdown   surface    is    expressed    in 

terms    of      the   set     of   polar      coordinates      R,   ij>,      $    (Figure      3).        For   a     point 

(R,    ip,    <}>)    on    the   breakdown'  surface,    using   Cartesian   coordinates    (X,    Y,      Z)    and 

the   angle   ui  between  X-axis    and      the    line-of-sight  Lgs  we   obtain      the   following 

geometrical   relationships: 

X   =    R  COSi^      ;      Y   =    (A0  +   R  Sllty)    COS<(>      ;      Z   =    (AQ  +   R  SINi|;)SIN<j>  ■ 

Ao 
C0S*MAX(M)    =      Cao+R(m)SIN^J  (3.3) 

U   =   U(C0S9,    SIN9    COS(|»,    SIN9    SIN<(>) 

Lbs/I^BsI  =  (COSoj  ,  SINo)  COSB,  SINu  SINS) 


The  cosines  COSk,  COSag,  COS<*s  in  (3.2;  are  expressed  as  scalar  products 
of  [LrS^ILBsN  an<*  unit  vectors  along  the  local  velocity  vector  U,  along 
the  local  normal  to  the  breakdown  surface  and  along  the  local  normal  to  the 
spacecraft  surface,  correspondingly. 

The  integration  is  performed  numerically  in  two  phases,  the  first  being 
the  integration  along  the  supplimentary  breakdown  surface  (Figure  4).  For 
this  first  phase,  the  straight  line  segment  which  constitutes  the 
supplementary  breakdown  surface  is  divided  into  several  intervals  of  length  AR 
(typically  10  intervals).  Each  interval  generates  a  half-strip  by  rotating  it 
from  <(>  =  0  through  $  =  ^(M  ).  This  strip  is  in  turn  subdivided  into 
several  sub- intervals  of  A<j>  each  (typically  10  intervals).  The  total  flux 
arriving  at  Xs  is  obtained  by  summing  contributions  from  each  sub-interval 
(two-dimensional  integration).  When  the  integration  along  the  supplementary 
breakdown  surface  is  concluded,  it  is  continued  into  the  breakdown  surface, 
where  AR  intervals  are  replaced  by  breakdown  surface  intervals  that  correspond 
to  a  fixed  Mach  number  increment  AM  (typically  AM=0.1).  The  integration 
proceeds  along  the  breakdown  surface  (2.5)  until  the  contribution  of  the  last 
AM  strip  is  negligibly  small.  The  computation  time  is  modest  (about  1  second 
CPU  per  Xs  point,  on  IBM  3033  mainframe  computer).  The  computations  were 
carried  out  by  a  code  RINGBD  written  specifically  for  this  purpose.  Further 
details  of  the  scheme  and  programming  can  be  obtained  by  reading  this  code 
which  is  given  in  Appendix  A. 


4.   RESULTS  AND  DISCUSSION 

4.1   Presentation  of  Results 

The  molecular  flux  backscattered  to  the  spacecraft  from  the  surface 
of  continuum  flow  breakdown  in  the  lip-centered  rarefaction  fan,  has  been 
computed  for  all  five  species  H,  HF,  H2  >  DF,  He.  The  results  are  depicted  in 
Figures  5  to  9  respectively.  For  each  species  two  more  cases  were  computed  in 
addition  to  the  nominal  case  (1.1),  where  the  stagnation  density  p0  was 
replaced  by  po/10  anc*  by  Po*^  (see  Figures  5  to  9).  This  has  been  done  in 
order  to  demonstrate  the  effect  of  variations  in  exit  flow  conditions  on  the 
flux.  The  particular  choice  of  p0  was  motivated  by  the  fact  that  the  effects 
of  changing  p0  are  not  obvious.  The  effects  of  changing  the  exit  Mach  number 
M  or  the  stagnation  temperature  T0  are  rather  obvious  (a  higher  flux  would 
result  from  either  a  decrease  in  M  or  an  increase  in  T0) .  It  turns  out  that 
for  points  lying  not  too  near  the  nozzle  lip  (Xs  >  0.1  m) ,  the  lower  density 
flow  generates  a  higher  backscattered  flux! 

In  addition  to  varying  pQ,  we  also  varied  the  breakdown  parameter  B, 
obtaining  a  surprising  result.  The  computation  was  performed  for  a  particular 
species  (HF) ,  and  the  results  obtained  upon  replacing  B=0.05  (nominal  value) 
by  B/2  and  by  B*2  are  brought  in  Figure  10. 

It  turns  out  that  the  B/2  case  has  the  higher  flux.  This  is  somewhat 
surprising,  since  a  lower  value  of  B  in  a  centered  rarefaction  fan  (equation 
2.5)  means  that  the  breakdown  of  continuum  flow  takes  place  in  a  region 
further  out  from  the  corner.  In  a  source  flow  (e.g.,  a  spherical  source), 
that  implies  lower  density  and  temperature,  which  would  give  rise  to 
lower  thermally  backscattered  flux. 
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An  explanation  to  these  seemingly  counterintuitive  results,  along 
with  some  deeper  insight  into  the  breakdown  surface  model  as  it  is  applied  to 
a  centered  rarefaction  flow,  can  be  obtained  by  taking  a  close  look  at  the 
flow  field  and  the  breakdown  surface  in  the  vicinity  of  the  corner.  We  take 
up  this  matter  in  the  following  sections. 

We  conclude  the  presentation  of  results,  by  comparing  the  flux  (in 
the  nominal  case)  of  the  five  species  with  each  other  (Figure  11).  This 
figure  underlines  the  fact  that  the  flux  of  light  species  (H,  H2 »  He)  is  many 
orders  of  magnitude  (typically  1015)  times  that  of  heavy  species  (HF,  DF) . 
Indeed,  these  results  demonstrate  a  well  known  effect:  When  an  expanding 
gaseous  mixture  of  light  and  heavy  molecules  experiences  a  breakdown  of 
continuum  flow,  a  separation  of  species  takes  place  (see  e.g.,  the  work  of 
Cattolica  et.  al.  [8]). 

4.2   The  Breakdown  Surface  and  Streamlines 

Consider  the  parametric  description  Rg(M)  for  the  breakdown  surface 
(Equation  2.5).   Normalizing  R  relative  to  the  exit  mean  free  path  \\ ,  we  get: 


RB(M)    =    RB(M1)[(M2    -    l)/(Ml    -    1)]  [(1    +JCiM2)/(l    +rl  Mj]" 


RgCM^/Xj   =    [(Ytt/2)1/2/(y+1)b]    (M*  -    1)1/2 

(4.1) 
1 

X^V'^Nof1    [1+^M2]    Y-l 


The  normalized  surface   Rg(M)/X   is   thus  independent   of   stagnation  density 
depending  only  on  y,  M  ,  B. 
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Let  us  now  derive  a  parametric  equation  Rg(M)  for  a  streamline  that 
enters  the  fan  at  point  Rg(M  )  on  the  exit  characteristic.  The  following 
geometrical  relationship  is  readily  obtained  by  considering  two  characteristic 

lines  \p  and  t|H-A^  and  a  streamline  inclined  at  the  Mach  angle  \i   to  them: 

dR«(\|/) 

— f =  -  Rg(i|>)  (tanu)-l  (4.2) 

ay 

Using   the  standard   Prandtl-Meyer   functions   (2.1),   we   get   the 

following  differential  equation  for   Rg(M): 

This  equation  is  readily  integrated,  giving: 


RS(M)  =  RgCl^)   [(1  +  I^-M2)/(l  +  £***)]    2(T~°  (4.4) 

As  pointed  out  by  Birdl^J,  there  is  a  particular  streamline  Rga(M) 
which  asymptotically  approaches  the  breakdown  surface  for  large  M,  since  the 
ratio  Rg(M)/Rg(M)  tends  to  a  constant  (not  zero)  when  M  +  ».  (Strictly 
speaking,  this  holds  only  for  hard-sphere  molecules,  i.e.,  only  when  cj  ■  0.5 
inf7b.      The   limit    is: 

lim  Rs(M>        W    r    m[-1_ //2 


RB(M)    =   K^("i)    I    2  2    1 

»  D  Mi    +  - 


Y" 

For  the  limiting  streamline  Rga(M) ,  the  ratio  Rga(M)/RB(M)  should 
tend  to  1.  This  determines  the  point  Rsa(M  )  at  which  the  limiting  streamline 
enters    the   fan,    as   well   as    the   entire   line    Rga(M) : 


Rsa(M1)    =    RB(Ml)       [    l      X1] 

M«  -  1 


2  2 
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Clearly,  Rsa(M)  is  larger  Chan  Rg(M)  for  any  14  >  M  ,  so  that  no 
streamline  beyond  Rsa(M)  can  cross  the  breakdown  surface.  This  pattern  is 
shown  in  Figure  4,  where  Rga(M)  is  denoted  "streamline  2",  and  "streamline  1" 
is  the  streamline  RgCM^  =  Rb(M,). 

All  this  leads  to  the  following  observation  regarding  the  continuum 
breakdown  of  the  flow  in  a  centered  rarefaction  fan^'K  Referring  to  Figure 
4,  the  fluid  entering  the  fan  through  the  supplementary  breakdown  surface 
(i.e.,  through  the  exit  characteristic  between  the  corner  and  the  initial 
point  of  streamline  1),  experiences  breakdown  immediately  upon  crossing  this 
surface.  Every  streamline  between  streamline  1  and  streamline  2  crosses  the 
breakdown  surface  at  some  Mach  number  M  >  M.  ,  and  at  that  point  the 
continuum  flow  breaks  down.  All  fluid  entering  the  fan  beyond  stream -line  2 
will  never  pass  through  the  breakdown  surface,  and  hence  will  maintain  a 
continuum  flow  regime  all  the  way  to  infinity.  Of  course,  that  is  only  true 
for  planar  centered  rarefaction  fans.  When  the  exhaust  flow  emerges  from  a 
nozzle  of  finite  width,  and  especially  when  the  exhaust  jet  has  a  ring 
symmetry  (as  in  our  case),  the  breakdown  surface  gradually  curves  in  a 
balloonlike  shape  towards  the  opposite  nozzle  lip,  forming  the  familiar  plume 
pattern  (Figure  1). 

4.3  Analysis  and  Discussion  of  Results 

The  foregoing  analysis  is  now  used  to  explain  the  variation  in  back- 
scattered  flux  due  to  a  change  in  exhaust  flow  conditions  at  the  nozzle  exit. 
Specifically,  we  consider  a  tenfold  decrease  in  stagnation  density'  (i.e.,  the 
case  po/10)>  an<3  hence  a  tenfold  increase  in  the  exit  mean  free  path  X  . 

The  effusion  flux  from  the  breakdown  surface  is  proportional  to  the 
local  density,  so  one  would  expect  to  observe  a  decrease  in  flux,   rather  than 
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an  increase  (see  Figures  5  to  9 ,  for  Xs  >  0.1  m) .  Other  factors  causing 
increased  flux,  must  then  be  larger  than  10  so  that  they  more  than  offset  the 
1/10  factor  in  density.  It  turns  out  that  these  effects  are  mainly 
geometrical,  in  that  a  tenfold  increase  in  \  causes  the  domain  of  integration 
on  the  breakdown  surface  to  increase  more  than  tenfold.  In  the  (X,Y)  plane 
there  is  a  tenfold  "blowup"  of  the  breakdown  surface,  due  to  the  self-similar 
structure  of  the  Prandtl-Meyer  flow  field.  As  a  result  of  this  "blowup"  in 
(X,Y),  the  angular  integration  range  <J>max  also  increases,  albeit  not  linearly 
(Equation  3.3).  Another  geometrical  effect  is  an  increase  in  the  flux 
incidence  cosine  factor  cosas  (see  Equation  3.2),  which  for  points  Xg 
sufficiently  far  from  the  nozzle  lip,  increases  roughly  tenfold  (while  the 
other  cosine  factor  cosag  is  almost  constant).  All  this  provides  a" 
qualitative  explanation  for  the  observed  increase  in  flux  at  far  points  (Xs  > 
0.1  m). 

As  for  the  near  range  (Xs  <  0.1  m) ,  another  effect  becomes 
increasingly  significant  as  Xs  approaches  the  nozzle  lip.  The  turning  angle 
k,  by  which  backscattered  molecules  have  to  be  deflected  relative  to  the  flow 
velocity  vector  in  order  to  reach  point  Xs  on  the  spacecraft  (Figure  2), 
increases  with  the  size  of  the  breakdown  surface  (fixed  M  and  Xs).  Since  the 
local  effusion  flux  (Equation  2.6)  decreases  rather  sharply  as  <  is  increased, 
the  net  result  is  a  tendency  to  get  a  reduced  backscattered  flux  at  near 
points  such  as  Xg  =  0.01  m  (Figures  5  to  9). 

We  now  turn  to  the  effect  of  changing  the  value  of  the  breakdown 
parameter  B.    From  equation  4.1  it  is  clear  that  multiplying  B  by  some  factor 
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will  have  the  same  "blowup"  effect  as  dividing  M  by  the  same  factor.  A 
tenfold  decrease  in  B  is  thus  geometrically  equivalent  to  a  tenfold  decrease 
in  pQ.  However,  since  the  local  effusion  flux  at  the  breakdown  surface  is 
proportional  to  p0  while  it  is  independent  of  B,  the  B/10  case  will  have  ten 
times  as  much  backscattered  flux  as  the  po/10  case.  In  order  to  illustrate 
the  sensitivity  of  the  flux  estimates  to  an  uncertainty  in  the  appropriate 
value  of  B,  we  computed  the  cases  B/2  and  B*2  for  one  species  (HF),  and  they 
are  presented  in  Figure  10.  The  variation  in  flux  relative  to  the  nominal 
case  (B  =  0.05),  is  by  a  factor  no  larger  than  about  5.  Results  for  other 
species  were  found  to  exhibit  comparable  variations. 

Does  this  observation  about  the  dependence  of  the   breakdown  surface 
on  B  agree  with   the  breakdown  surface   appropriate  to   the  far  field   of  the 

exhaust  plume?    In  stationary  source  flow  into  vacuum,  and  when  M  »  1,  the 

5-1 
breakdown  parameter  varies   with  radius  as   B  ~  R    (5=1  for  cylindrical 

source,  5=2  for  spherical  source).   In  a  ringjet,  the  stream  tubes  of  the 

exhaust  plume  generally  diverge  at  a  rate  higher  than  that  of  stream   tubes  in 

a  cylindrical  source   flow,  so  the  effective  value  of  5  in  a  ringjet  is  5  >  1 . 

Hence,  in   this  case  the   far  field   breakdown  surface  moves   downstream  along 

each  stream  tube  as   the  value  of  B  increases.    This  is   indeed  geometrically 

compatible  with  the  fact  that  near  the  corner  of  the   lip-centered  rarefaction 

_1 
fan  B  ~  R   ,  as  shown  schematically  in  Figure  12.  The  dependence  of  the  break- 
down surface  on  B  near   the  corner  and   in  the   far  field,  thus   assures  that 
complete   breakdown  surfaces  corresponding   to  different  values  of  B,   do  not 
intersect  (Figure  12). 

In   the  foregoing  discussion   it  was   pointed  out  that  variations  in 
flux  caused  by  changes  in  parameters  such  as   p0  and  B,  were   directly  related 
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to  the  self-similar  structure  of  the  Prandtl -Meyer  flow  field.  It  has  been 
further  shown  that  these  variations  are  well-understood  within  the  framework 
of  the  breakdown  surface  model  and  that  they  are  not  excessively  large.  Are 
we  to  conclude  that  the  thermally  backscattered  flux  estimates  of  the  present 
model  are  also  physically  plausible  and  reliable?  In  the  following  section  we 
take  up  this  matter,  arriving  at  some  interesting  conclusions  about  this  model 
and  its  range  of  validity. 
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4.4  C r itical  Examination  of  the  Model 

Consider  the  centered  rarefaction  flow  field  of  a  compressible  fluid 
negotiating  an  expansive  corner  at  supersonic  speed  (Prandtl-Meyer  flow).  The 
streamlines  of  this  flow  field  have  an  orderly  "layered"  structure,  with  each 
streamline  curving  around  the  corner,  starting  at  its  point  of  entrance  into 
the  fan  (see  Figure  4). 

The  present  model  is  based  on  the  stipulation  that  there  is  a  point 
of  continuum  flow  breakdown  on  each  streamline,  provided  this  streamline  is 
not  beyond  a  certain  limiting  streamline.  Consider  a  sample  molecule  effusing 
from  this  breakdown  point  toward  the  spacecraft.  It  advances  at  constant 
speed  along  a  straight  line  trajectory,  traversing  all  inner  streamlines. 
Since  the  flow  velocity  vector  points  away  from  the  spacecraft,  and  since  the 
flow  is  highly  supersonic  so  that  the  velocity  of  most  individual  molecules 
does  not  differ  much  from  the  flow  velocity  (i.e.,  it  is  a  "cold"  flow),  any 
collision  of  the  sample  molecule  with  a  mainflow  molecule  will  most  probably 
divert  the  sample  molecule  away  from  the  spacecraft..  What  is  the  probability 
that  a  sample  molecule  would  traverse  this  cross  flow  collisionlessly?  This 
probability  is  simply  exp(-n),  where  n  is  the  expected  number  of  collisions 
along  the  straight-line  trajectory  from  the  point  of  breakdown  to  the 
spacecraft.  In  the  typical  operating  conditions  assumed  here,  we  estimated  n 
to  be  roughly  about  10.  Since  this  no-collision  probability  factor  is  ignored 
in  the  formulation  of  the  present  model,  the  backscattered  flux  may  be 
exaggerated  by  a  factor  of  exp(10)  or  about  101*.  We  conclude  that  in  all 
likelihood,  the  prediction  of  the  breakdown  surface  model  for  thermally 
backscattered  flux  from  a  centered  rarefaction  flow,  is  substantially 
overestimated. 
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Can  anything  be  done  to  improve  the  present  model?  One  may  be 
inclined  to  suggest  at  this  point  that  the  obvious  remedy  is  to  incorporate 
the  no-collision  probability  factor  into  the  model.  Rather,  we  prefer  to 
retain  the  breakdown  surface  model  in  its  present  form  as  a  simple  means  of 
obtaining  an  overestimate  to  the  thermally  backscattered  flux  from  a  centered 
rarefaction  flow.  An  improved  model  can  be  constructed  by  considering  thermal 
backscattering  from  the  entire  flow  field  (tempered  by  the  probability  of 
no-collision),  without  resorting  to  the  physically  untenable  notion  of  an 
abrupt  transition  from  continuum  flow  to  free  molecular  flow. 
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5.   CONCLUSIONS 

Some  surprising  similarity  laws  of  the  breakdown  surface  model  were 
observed.  It  has  been  shown  that  they  were  a  direct  result  of  the  self 
similar  structure  of  the  Prandtl-Meyer  flow  field  to  which  the  model  was 
applied.  Specifically,  it  was  found  (and  shown  plausible)  that  reduced  values 
of  either  the  exhaust  stagnation  density  p0»  or  the  breakdown  parameter  B, 
caused  higher  backscattered  flux. 

The  breakdown  surface  model  for  thermally  backscattered  flux  from  a 
centered  rarefaction  fan,  has  been  shown  to  overestimate  the  flux  arriving  at 
the  spacecraft.  It  is  suggested  that  an  improved  model  be  constructed  by 
considering  thermal  backscattering  from  the  entire  flow  field,  along  with  the 
probability  factor  for  a  side-scattered  molecule  traversing  the  main  flow 
collisionlessly. 

The  molecular  flux  of  corrosive  species  (HF,  DF)  arriving  at  the 
spacecraft  (Figures  6  and  8)  is  no  larger  than  about  107  (sec-  m~2),  which  is 
negligible  since  it  corresponds  to  about  10  molecular  monolayers  per  year. 
This  conclusion  is  reliable  since  even  this  flux  level  is  an  overestimate. 

The  maximum  thermally  backscattered  flux  of  light  species  (H,  H2,  He)  is 
in  the  range  of  1020  to  1022  (sec-1  m"2)  (see  Figures  5,  7,  9).  Thus,  we 
conclude  that  while  thermal  backscattering  would  contribute  significantly  to 
the  flux  of  light  molecules  arriving  at  the  spacecraft,  it  is  utterly 
negligible  as  far  as  heavy  molecules  are  concerned. 
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APPENDIX  A.   The  Computer  Code  RINGBD 

We  present  a  printout  of  the  code  RINGBD  along  with  the  reults  of  the 
nominal  case  (printout  of  actual  run).  This  is  preceded  by  a  brief 
description  of  the  subroutines  and  a  summary  of  major  variables  with  their 
code  and  report  notations. 

A.l   Description  of  Subroutines 


MAIN  PROGRAM 


INIDAT 


FLUX 


BREAKR 


BREAKM 


Computes  flux  integration  by  summation  of  segment 
contributions  (centered).   Printing  of  results. 

-  Definition  of  all  data  (no  input  file). 
Preparatory  evaluation  of  parameters.  Printing 
of  data. 

Evaluates  flux  emitted  from  a  single  point  on 
breakdown  surface  (mean  segment  values)  to  a 
point  on  spacecraft  (XS). 

Computes  point  on  breakdown  surface  for  given 
Mach  number. 

-  Computes  point  on  breakdown  surface  for  mean 
Mach  number  of  a  segment. 
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A. 2   Code  Versus  Report  Notation 


XC(I)  -  [AB]  - 

WC(I)  -  Wi 

WAV   -  WA 

TO    -  T0 

RHOO  -  p0 

G     -  Y 

EMI    -  M, 


LAMDA1-  X 


AO     -  Aq 
R     -  R 
DIST  -  LBS 


xs 

"   Xg 

PSI 

-  ,p 

AMU 

-  U 

TETA 

-  9 

PHI 

"  4> 

W 

-  O) 

BETA 

-  6 

DMO    -  AM 
EM    -  M 
PBIRD  -  B 


Mole  fraction  of  species  AB.    (1=1,2,3,4,5  cor- 
responds to  H,  HF,  H2,  DF,  He). 
Molecular  weight  of  species  1. 
Average  molecular  weight 
Stagnation  temperature 
Stagnation  density 
Specific  heat  ratio 
Exit  Mach  number 
Exit  mean  free  path 
Spacecraft  radius 

Distance  from  corner  (X=0,  Y2+Z2  =  Aq2). 
Distance  between  emitting  point  on  breakdown  surface 
and  receiving  point  (XS)  on  spacecraft. 
Point  on  spacecraft  (X=Xg,  Y=AQ,  Z=0). 
Characteristic  angle 
Mach  angle 

Velocity  vector  angle 
Rotation  angle  for  flux  integration 

Angle  between  x-axis  and  line-of -sight  Lyg 

-♦■ 
Angle   between  Y-axis  and  projection  of  Lgs  on  (Y,Z) 

plane. 

Mach  number  increment  for  flux  integration 

Mach  number 

Breakdown  parameter 
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A.  3  Code  Listing  (Run  of  Nominal  Case) 


$J0B 


6 

7 

8 

9 

10 

11 


12 


101 


RINGBD,NOXREF 

IMPLICIT  REAL*3(A-H,0-Z,$) 

COMMON  /GAMA/G,G1,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14 
1  G16,G17,G18,G19,G20 

COMMON  /PAR/CO, ENO, EMI, D, TL IM, ETAL IM, CLIM, EL 0,Q0, TO, 
1  PBIRD,RBIRD, DM0, DEG, OMEGA, XSV( 51) 

COMMON  /NPAR/NETA,NC,NT,NEM0,NPHI,NXS,NR0,NSPEC 

COMMON  /GE0M/APF,PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1, 

1  CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF, 

2  PHI,SPHI,CPHI,RMIN,RMAX,XS,DIST, 

3  AMU1,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0 
COMMON  /EPSIL/EPSQ,EPSETA,EPST,EPSC,EPSEM 

COMMON  /EXTREM/TEXT , ETAEXT , CEXT , REXT , PSI EXT , EMEXT , BEXT , QEXT 
COMMON  /SPEC/WAV,XC(5),WC(5),WCR(5),XNAME(5),QC(5),FLUXC(5) 
DIMENSION  DSUMC5) 
PRINT  101 

FORMATCl'/lX, 'RINGBD  -  FLUX  INTEGRATION  FROM  BREAKDOWN', 
1  IX, 'SURFACE'//) 


CALL  INIDAT 


13 

PRINT  110,XNAME 

14 

110 

FORMATC///1X, '   NX','  NEM','     XS    ', 
1            5(4X,A6,1X, '/  LOG', IX)) 

15 

DO  200  NX=1,NXS 

16 

EM=EM1 

17 

CALL  BREAKR(EM,RF) 

18 

IF(NRO.GT.O)  RF=R0 

19 

XF=RF*CPSI1 

20 

YF=RF*SPSI1+A0 

21 

XS=XSV(NX) 

22 

QMAX=0. 

23 

DO  45  N=1,NSPEC 

24 

FLUXC(N)=0. 

25 

45 

CONTINUE 

26 

DO  1  NEM=1,NEM0 

27 

RN  =  RF 

28 

XN  =  XF 

29 

YN=YF 

30 

IF(NEM.GT.NRO)  GO  TO  41 

31 

RF=R0+DFLOAT(NEM)X(RMIN-R0)/DFLOAT(NR0) 

32 

XF=RFXCPSI1 

33 

YF=RF*SPSI1+A0 

34 

RMEAN=(RN+RF)/2.D0 

35 

EMMEAN=EM1 

36 

PSIM=PSI1 

37 

SPSIM=SPSI1 

38 

CPSIM=CPSI1 

39 

GO  TO  42 

40 

41 

CONTINUE 

41 

EM=EM+DM0 

42 

CALL  BREAKR(EM,RF) 

43 

EMMEAN=EM-DM0/2.D0 

44 

CALL  BREAKM(EMMEAN,RMEAN) 

45 

42 

CONTINUE 

46 

AL0NG=DSQRT((XF-XN)XX2+(YF-YN)XX2) 

47 

SALFA=(YF-YN)/ALONG 

48 

CALFA=(XF-XN)/ALONG 

49 

PHIMAX=DARCOS(A0/(A0+RMEAN*SPSIM)) 

50 

DPHI=PHIMAX/NPHI 

PHIMAX' 


QMAX 


RI 
RI 
,G15,RI 
RI 
■  RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 
RI 


N00010 
N00020 
N00030 
N00040 
N00050 
N00060 
N00070 
N00080 
N00090 
N00100 
N00110 
N00120 
N00130 
N00140 
N  0  0 1 5  0 
N00160 
N  0  0 1 7  0 
N00180 
N00190 
N00200 
N00210 
N00220 
N00230 
N00240 
N00250 
N00260 
N00270 
N00280 
N00290 
NO  030  0 
N00310 
N00320 
N00330 
N00340 
N00350 
N00360 
N00370 
N00380 
N00390 
N00400 
N00410 
N00420 
N00430 
N00440 
N00450 
N00460 
N00470 
N00480 
N00490 
N00500 
N00510 
N00520 
N00530 
N00540 
N00550 
N00560 
N00570 
N00580 
N00590 
N00600 
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51  DO  44  N  =  1,NSPEC  RIN0061 

52  DSUM(N)=0.  RIN0062 

53  44  CONTINUE  RIN0063 

54  DO  2  NP=1,NPHI  RIN0064 

55  PHI=(DFL0AT(NP)-0.5D0)XDPHI  RIN0065 

56  CALL  FLUX(EMMEAN,RMEAN)  RIN0066 

57  CROSSl=SW*CBETA  RIN0067 

58  CR0SS2=(SALFA)x(-CW)+(-CALFAXCPHI)x(-SWXCBETA)+  RIN0068 
1        (-CALFA*SPHI)X(-SW*SBETA)  RIN0069 

59  GOREM=CROSSlxCROSS2xDPHlx(A0+RMEANXSPSIM)xALONG/DISTXX2  RIN0  07  0 

60  DO  24  N=1,NSPEC  RIN0071 

61  DSUM(N)=DSUM(N)+QC(N)*GOREM  RIN0072 

62  24  CONTINUE  RIN0073 

63  IF(QMAX.GE.QEXT)  GO  TO  25  RIN0074 

64  QMAX=QEXT                           .  RIN0075 

65  25  CONTINUE  RIN0076 
C  PRINT  22,NEM,NP,EMMEAN,RMEAN,PHIMAXXDEG,DARC0S(CW)XDEG,  RIN0077 
C  1  BETAXDEG,PHIXDEG,CR0SS1,CR0SS2, ALONG, DIST,GOREM,  RIN0078 
C  2  SPSIM,A0,(QC(N),FLUXC(N),N=1,NSPEC)  RIN0079 
C22  F0RMAT(/1X, « NEM, NP, EMMEAN, RMEAN, PHIMAX= ' , 213, 3D13 . 4/  RIN0080 
C  1  IX, 'W,BETA,PHI=',3D15.5/  RIN0081 
C  2  IX, •CR0SS1,CR0SS2,AL0NG,DIST,G0REM=,,5D15.5/  RIN0082 
C  3  IX, 'SPSIM,A0=',2D15.5/  RIN0083 
C  4        IX, ,QC,FLUXC=',5(1X,D10.3,1X,D10.3))  RIN0084 

66  2  CONTINUE  RIN0085 

67  DO  26  N=1,NSPEC  RIN0086 

68  FLUXC(N)=FLUXC(N)+DSUM(N)  RIN0087 

69  26  CONTINUE  RIN0088 

70  IF(NEM.LE.NR0+2)  GO  TO  1  RIN0089 

71  DO  27  N=1,NSPEC  RIN0090 

72  IF((DSUM(N)/FLUXC(N)).GT.EPSEM)  GO  TO  28  RIN0091 

73  27  CONTINUE  RIN0092 

74  GO  TO  10  RIN0093 

75  28  CONTINUE  RIN0094 

76  1  CONTINUE  RIN0095 

77  10  CONTINUE  RIN0096 

78  DO  31  N=1,NSPEC  RIN0097 

79  FLUXC(N)=2.D0XXC(N)*FLUXCCN)  RIN0098 

80  31  CONTINUE  RIN0099 

81  PRINT  11,NX,NEM,XS,PHIMAXXDEG,QMAX,  RIN0100 
1          (FLUXC(N),DABS(DLOG10(FLUXC(N))),N=1,NSPEC)  RIN0101 

82  11  F0RMAT(/1X,2I4,F9.4,F7.2,D11.3,5(1X,D10.3, V',F5.2))  RIN0102 

83  200  CONTINUE  RIN0103 

84  PRINT  102  RIN0104 

85  102  FORMATC///1X, 'END  RINGBD  RUN',///)  RIN0105 

86  STOP  RIN0106 

87  END  RIN0107 

88  SUBROUTINE  INIDAT  RIN0108 

89  IMPLICIT  REALX8(A-H,0-Z,$)  RIN0109 

90  REALX8  LAMDA0,LAMDA1  RIN0110 

91  COMMON  /GAMA/G,G1,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14,G15,RIN0111 
1             G16,G17,G18,G19,G20  RIN0112 

92  COMMON  /PAR/CO, ENO, EMI, D, TLIM, ETALIM, CLIM, ELO, QO , TO,  RIN0113 
1            PBIRD,RBIRD, DMO, DEG, OMEGA, XSV(51)  RIN0114 

93  COMMON  /NPAR/NETA,NC,NT,NEMO,NPHI,NXS,NRO,NSPEC  RIN0115 

94  COMMON  /GE0M/APF,PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1,  RIN0116 

1  CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF,  RIN0117 

2  PHI,SPHI,CPHI,RMIN,RMAX,XS,DIST,  RIN0118 

3  AMU1,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0  RIN0119 
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95  COMMON  /EPSIL/EPSQ,EPSETA,EPST,EPSC,EPSEM  RIN01200 

96  COMMON  /EXTREM/TEXT,ETAEXT,CEXT,REXT,PSIEXT,EMEXT,BEXT,QEXT  RIN01210 

97  COMMON  /SPEC/WAV, XC( 5) , WC( 5) , WCR( 5) , XNAMEC5) , QC(5) , FLUXC(5)  RIN01220 

98  DATA  XC/.091D0, .091D0, .104D0, .135D0, .579D0/  RIN01230 

99  DATA  WC/1.00D0,20.0D0,2.00D0,21.0D0,4.00D0/  RIN01240 

100  DATA  XNAME/'    H   ',»   HF   »,»   H2   ','   DF   »,«   HE   V  RIN01250 

101  PAI=4.XDATAN( .ID  1)  RIN01260 

102  AR=8.3143D3  RIN01270 

103  AV=6.022D  26  RIN01280 
C  0MEGA=0.5  IS  FOR  HARD  SPHERE  COLLISIONS,  RIN01290 
C  AN  AVERAGE  RECOMMENDED  VALUE  IS  ABOUT  OMEGA  =  0.75  RINO-1300 

10<^  0MEGA  =  0.5D0  RIN01310 

105  NSPEC=5  RIN01320 

106  WAV=0.  RIN01330 

107  DO  51  N=1,NSPEC  RIN01340 

108  WAV=WAV+XC(N)*WC(N)  .  RIN01350 

109  51    CONTINUE  RIN01360 

110  DO  52  N=1,NSPEC  RIN01370 

111  WCR(N)=DSQRT(WC(N)/WAV)  RIN01380 

112  52    CONTINUE  RIN01390 

113  A0=2.5D0  RINOl^OO 

114  EM1=4.0D0  RIN01410 

115  RHO0=0.0075D0  RIN01420 

116  T0=2.300D3  RIN01430 

117  D=2.50D-10  RIN01440 

118  G=1.54D0  RIN01450 

119  .  EN0=RHO0*AV/WAV  RIN01460 

120  CO=DSQRT(GXARXTO/WAV)  RIN01470 

121  PBIRD=0.05D0X2.D0  RIN01480 
C  RO  IS  THE  RADIUS  FOR  BEGINNING  THE  INTEGRATION  ALONG  THE  M=M1  RIN01490 
C  CHARACTERISTICCTHE  AUGMENTED  BREAKDOWN  SURFACE).  RIN01500 
C  NRO  IS  THE  NUMBER  OF  INTEGRATION  INTERVALS  ON  THIS  SEGMENT.  RIN01510 
C   FOR  NO  INTEGRATION  ALONG  M=M1  CHARACTERISTIC,  SET  NR0=0.  RIN01520 

122  R0=0.  RIN01530 

123  NR0=10  RIN01540 

124  DM0=0.1D0  RIN01550 

125  NEM0=20.D0/DM0+NR0  RIN01560 
C  TO  GET  FLUX  DUE  TO  AUGMENTED  BREAKDOWN  SURFACE  SOLELY,  ACTIVATE:  RIN01570 
C  NEM0=NR0  RIN01580 

126  NPHI=10  RIN01590 

127  NXS=13  RIN01600 

128  XSI=l.D-2  .  RIN01610 

129  XSF=1.D1  RIN01620 

130  XSV(1)=XSI  RIN01630 

131  IF(NXS.EQ.l)  GO  TO  111  RIN0164Q 

132  DXL=(DLOG(XSF)-DLOG(XSI))/(NXS-l.D0)  RIN01650 

133  XLI=DLOG(XSI)  RIN01660 

134  DO  11  NX=2,NXS  RIN01670 

135  XSV(NX)=DEXP(XLI+(NX-1.D0)XDXL)  RIN01680 

136  11    CONTINUE  RIN01690 

137  111   CONTINUE  RIN01700 

138  EPSEM=l.D-5  RIN01710 

139  DEG=180.D0/PAI  RIN01720 

140  PAI2=PAI/2.D0  RIN01730 

141  GAMMA=G  RIN01740 

142  G1=(G-1.D0)/2.D0  RIN01750 

143  G2=(G+1.D0)/(2.D0X(G-1.D0))  RIN01760 

144  G3=G/2.D0  RIN01770 

145  G4=(G+1.D0)/(G-1 .DO)  RIN01780 

146  G5=DSQRT((G+1.D0)/(G-1.D0))  RIN01790 
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147 
148 
149 

150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 

168 
169 

170 
171 
172 
173 
174 
175 


176 
177 


178 
179 

180 
181 


201 


202 
203 
21 


22 


23 


G6=1.D0/(G-1.D0) 
G7=2.D0/(G+1.D0) 

G8=(5.D-lx(G+l.D0)**2/(G-l.D0))**(l,D0/CG+l.D0))  * 
1    ((G+l.D0)/(G-l.D0))xx((G-l.D0)/(G+l.D0)) 
G9=(G+3.D0)/(2.D0X(G-1.D0)) 
G10=(7.DO-3.DOXG)/(2.DOX(G-1.DO)) 
G11=DSQRT(G/PAI)/(2.D0X(G+1.D0)) 
G12=DSQRT(G/2.D0) 
G13=1.D0/DSQRT(2.D0XGXPAIXX3) 
LAMDA0=1.D0/(DSQRT(2.D0)XPAIXDXX2XEN0) 
LAMDA1=LAMDAOX(1.DO+G1XEM1XX2)XX(G6-OMEGA+0.5DO) 
RBIRD=G11/(DXX2XEN0XPBIRD) 
ZETA1=G5XDATAN(DSQRT(EM1XX2-1.D0)/G5) 
AMU1=DARSIN(1 .D0/EM1) 
PSI1=PAI2+AMU1 
SPSI1=DSIN(PSI1) 
CPSI1=DC0S(PSI1) 
PSIF=PAI2+AMU1+ZETA1-G5XPAI2 
SPSIF=DSIN(PSIF) 
CPSIF=DCOS(PSIF) 
CALL  BREAKR(EM1,RMIN) 
RSMIN=RMINXDSQRT((2.D0/(G-1.D0)+EM1XX2)/(EM1XX2-1.D0)) 


PRINT  201,NSPEC,XNAME 
F0RMAT(/1X, 'SPECIES  DATA 
1         IX, 'SPECIES  NAMES 
PRINT  202, XC 

FORMATC  IX, 'MOLE  FRACTION 
PRINT  203, WC 
FORMAK  1X,'M0L.  WEIGHT 
PRINT  21,AR,AV,WAV,G,RH001 
F0RMAT(/1X, 'THERMODYNAMIC 


NSPEC=',I3/ 

•,11(2X,A6,2X)) 

XC=*,11(F8.4,2X)) 


WC=',11(F8.4,2X)) 

T0,EN0,C0,D 

DATA'/ 

1  IX, 'AR,AV,WAV,GAMMA=*,2X,2D14.5,2F9.3/ 

2  IX, 'RHO0,T0,EN0,C0,D=',D12.4,F8.0,D13.5,2D12.4) 
PRINT  22,EM1,PSI1XDEG,PSIFXDEG,PBIRD, 

1  A0,RMIN,RSMIN,R0, 

2  LAMDA0,LAMDA1,RMIN/LAMDA1,RSMIN/LAMDA1 
F0RMAT(/1X, 'FLOW  AND  GEOMETRY  DATA'/ 

1  IX, »EM1,PSI1,PSIF,PBIRD=',4F9.3/ 

2  IX, »A0,RMIN,RSMIN,R0=',F9.3,3D14.5/ 

3  IX, 'LAMDA0,LAMDA1,RMIN/LAMDA1,RSMIN/LAMDA1=',4D13 
PRINT  23,DM0,NPHI,NR0,EPSEM 

FORMATC/1X, 'INTEGRATION  DATA'/ 
1         IX, 'DM0,NPHI,NRO,EPSEM=«,F9.3,2I5,4X,D12.3) 
RETURN 
END 


4) 


RINOli 
RINOli 
RINOli 
RINOli 
RINOli 
RINOli 
RINOli 
RINOli 
RINOli 
RINOli 
RIN01' 
RIN011 
RIN011 
RIN011 
RIN011 
RIN011 
RIN011 
RIN011 
RIN01' 
RIN01" 
RIN02I 
RIN02I 
RIN02I 
RIN02I 
RIN02I 
RIN02I 
RIN02I 
RIN02I 
RIN02I 
RIN02I 
RIN02! 
RIN02. 
RIN02! 
RIN02: 
RIN02! 
RIN02! 
RIN02! 
RIN02! 
RIN02: 
RIN02: 
RIN02; 
RIN02; 
RIN02; 
RIN02; 
RIN02; 
RIN022 


182 
183 
184 

185 

186 
187 


188 
189 


SUBROUTINE    FLUX(EM,R) 

IMPLICIT    REALX8(A-H,0-Z,$) 

COMMON  /GAMA/G,G1,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14,G15, 
1  G16,G17,G18,G19,G20 

COMMON  / PAR/CO, ENO, EMi , D, TLIM, ETAL IM, CLIM, EL 0,Q0, TO, 
1  PBI RD, RBI RD, DMO, DEG, OMEGA, XSV( 51) 

COMMON  /NPAR/NETA,NC,NT,NEMO,NPHI,NXS,NRO,NSPEC 

COMMON  /GE0M/APF,PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1, 

1  CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF, 

2  PHI,SPHI,CPHI,RMIN,RMAX,XS,DIST, 

3  AMU1,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0 
COMMON  /EXTREM/TEXT,ETAEXT,CEXT,REXT,PSIEXT,EMEXT,BEXT,QEXT 
COMMON  /SPEC/WAV,XC(5),WC(5),WCRC5),XNAME(5),QC(5),FLUXC(5) 


RIN02; 
RIN02J 
RIN022 
RIN022 
RIN02; 
RIN02] 
RIN02: 
RIN02: 
RIN02: 
RIN02; 
RIN02: 
RIN02: 
RIN02: 
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190  SPHI=DSINCPHI)  RIN02390 

191  CPHI=DCOS(PHI)  RIN02400 

192  XMEAN=RXCPSIM  RIN02410 

193  YMEAN=(A0+RXSPSIM)XCPHI  RIN02420 

194  ZMEAN=(AO+RXSPSIM)xSPHI  RIN02430 

195  TBETA=ZMEAN/(YMEAN-AO)  RIN02440 

196  BETA=PAI2-DATAN(1.D0/TBETA)  RIN02450 

197  SBETA=DSIN(BETA)  RIN02460 

198  CBETA=DCOS(BETA)  RIN02470 

199  DIST=DSQRT((XS-XMEAN)XX2+(YMEAN-A0)XX2+ZMEANXX2)  RIN02480 

200  CW=-(XS-XMEAN)/DIST  RIN02490 

201  SW=DSQRT(1 .D0-CWXX2)  RIN02500 

202  GM=(1 .D0+GlxEMXX2)xx(-G2)  RIN02510 

203  AMU=DARSINC1 .DO/EM)  RIN02520 

204  TETA=PSIM-AMU  RIN02530 

205  STETA=DSIN(TETA)  RIN02540 

206  CTETA=DCOS(TETA)  RIN02550 

207  CKAPA=(CTETA)X(-CW)+(STETAxCPHI)x(-SWXCBETA)+  RIN02560 
1       CSTETA*SPHI)X(-SWXSBETA)  RIN02570 

208  SKAPA=DSQRT(1.D0-CKAPAXX2)  RIN02530 

209  QEXT=0.  RIN02590 

210  DO  1  N=1,NSPEC  RIN02600 

211  EMT=EMXCKAPAXG12XWCR(N)  RIN02610 

212  IF(DABSCEMT) .GT.13.D0)  EMT=EMTX( 13 . DO/DABSC EMT) )  RIN02620 

213  P0N=G3XEMXX2XWCR(N)XX2  RIN02630 

214  P0WT=P0WXSKAPAX*2  RIN02640 

215  IF(POW.GT.1.0D2)POW=1.0D2  RIN02650 

216  IFCPOWT.GT.1.0D2)  POWT=1.0D2  RIN02660 

217  EXP1=DEXP(-P0W)  RIN02670 

218  EXP2=DEXP(-P0WT)  RIN02680 

219  ERFC1=DERFC(-EMT)  RIN02690 

220  IFCERFCl.LT. l.D-43)  ERFCl=l.D-43  RIN02700 
C  IF(XS.LT.5.D0)  GO  TO  234  RIN02710 
C  IF(EXPl.GT.l.D-20.AND.ERFCl.GT.l.D-20)  GOTO  234  RIN02720 
C  PRINT  235,N,P0W,P0WT,EMT,EXP1,EXP2,ERFC1  RIN02730 
C235  F0RMAT(I4,8D11.3)  RIN02740 
C234  CONTINUE  RIN02750 

221  EVER1=EN0X(C0/WCR(N))XG13XGMX(1.D0+EMTXX2)XEXP1  RIN02760 

222  EVER2=(EN0X(C0/WCR(N))X0.5D0/PAI)X(EMXWCR(N))XCKAPAXGMX  RIN027  7  0 
1        (1.5D0+EMTXX2)XEXP2XERFC1  RIN02780 

223  QC(N)=EVER1+EVER2  RIN02790 

224  IF(QEXT.GE.QCCN))  GO  TO  1  RIN02800 

225  QEXT=QC(N)  RIN02810 

226  1  CONTINUE  RIN02820 

227  RETURN  RIN02830 

228  END  RIN02840 

229  SUBROUTINE  BREAKRC EM, R)  RIN02850 

230  IMPLICIT  REALX8(A-H,0-Z,$)  RIN02860 

231  COMMON  /GE0M/APF,PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1,  RIN02870 

1  CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF,  RIN02880 

2  PHI,SPHi,CPHI,RMIN,RMAX,XS,DIST,  RIN02890 

3  AMU1,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0  RIN02900 

232  COMMON  /PAR/CO, ENO , EMI , D, TL IM, ETALIM,CLIM, EL  0, QO , TO ,  RIN02910 
1            PBIRD,RBIRD, DMO, DEG, OMEGA, XSVC51)  RIN02920 

233  COMMON  /GAMA/G, Gl ,G2, G3, G4, G5, G6 , G7 , G8 , G9, G10, Gil , G12, G13,G14,G15, RIN029  3  0 
1             G16,G17,G18,G19,G20  RIN02940 

234  R=RBIRDXDSQRT(EMXX2-1 .DO)x(l . D0+GlxEMXX2)xx(G6-OMEGA+0 . 5D0 )  RIN02950 

235  ZETA=G5XDATAN(DSQRT(EMX*2-1.D0)/G5)  RIN02960 

236  PSI=PAI2+AMU1+ZETA1-ZETA  RIN02970 
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237  XF=R*DCOSCPSI)  rt!!o2990I 

238  YF=R.DSIN(PSI)+AO  rJnoIoOo! 

HI            1           CONTINUE  RIjjjJIJlJ. 

2«  RETURN  .                                                           RIM03030 

242  ENTRY   BREAKMCEM,R)  RIN03050 

243  r=rbIRD*DSQRT(EM**2-1.D0)*C1.D0+G1*EMXX2)x*CG6-OMEGA+0.5D0)  g ? M 2 1 S f  2 

244  ZETA  =  G5*DATANCDSQRTCEM*.*2-1.D0)/G5)  5tSS?S5S 

245  PSIM=PAI2+AMU1+ZETA1-ZETA  rtKo3090 

246  SPSIM=DSIN(PSIM)  RIN03100 

247  CPSIM=DCOS(PSIM)  RIN03  10 

248  RETURN  RIN03120 

249  END  RINOil^U 


$ENTRY 


RIN03130 


HF 

H2 

DF 

ME 

0.0910 

0.1040 

0.1350 

0.5790 

20.0000 

2.0000 

21.0000 

4  0000 

■UNOBO  -  FLUX  INTEGRATION  FROM  BREAKDOWN  SURFACE 

SPECIES  DATA    NSP£C=   5 
SPECIES  NAMES  H 

MOLE  FRACTION  XC-       0.0910 
MOL.  HEIGHT    HC=   1.0000 

THERMODYNAMIC  DATA 

AR.AV,WAV,GAMMA=      0.83143D  04    0.60220D  27     7.270     1 . 540 

RHO0.T0.ENO,C0.D=   0.7500D-02    2300.   0.62125D24   0.2013D04   0.2500D-09 

FLOW  AND  GEOMETRY  DATA 

EM1.PSI1.PSIF.PBIRD=     4.000   104.478    41.044     0.100 

A0.RMIN,RSM1N,R0=     2.500    0  30374D-02    0.34812D-02    0.00000D00 

tAMOAO,LAMDAl,RMIN/LAMDAl,RSMIN/LAMDAl=    0.5797D-05    0.12810-03    0.2J72D  02    0.2718D  02 

INTEGRATION  DATA 

DMO.NPHI.NRO.EPSEM*     0.100    10    10        0. 1000-04 

XS  PHIMAX  flMAX           H    /  LOG        HF    /  LOG        H2    /  LOG        DF    /  LOG        HE    /  LOG 

0.0100  11.19  0.120D  25  0 .8240  22/21 . 92  0 .2690  07/  6 . 43  0.777D  21/20.89  0.6910  06/5.84  O.653D20/19.81 

0.0178  12.15  0.1210  25  0.474D  22/21.68  0.1620  07/  6.21  0.4240  21/20.63  0.4180  06/  5.62  0.3490  20/19.54 

0.0316  13.50  0.1220  25  0.237D  22/21.38  0.7760  06/  5.89  0.1960  21/20.29  0.202D  06/  5.30  0.154D  20/19.19 

0.0562  15.27  0.123D  25  0.1060  22/21.02  0.2950  06/  5.47  0.780D  20/19.89  0.7700  05/  4  89  0.5680  19/18.75 

0.1000  17.16  0.1240  25  0.4210  21/20.62  0.8920  05/  4.95  0.268D  20/19.43  0.2330  05/  4.37  0.175D  19/18.24 

0.1778  19.55  0.124D  25  0.1490  21/20.17  0.219D  05/  4.34  0.791D  19/18.90  0.575D  04/  3.76  0.451D  18/17.65 

0.3162  22.52  0.1240  25  0.471D  20/19.67  0.464D  04/  3.67  0.203D  19/18.31  0.1220  04/  3.09  0.100D  18/17.00 

0.5623  25.65  0.124D  25  0.1320  20/19.12  0.9010  03/  2.95  0.4610  18/17.66  0.2370  03/  2.37  0.2010  17/16.30 

1.0000  29.39  0.124D  25  0.3350  19/18.52  0.1670  03/  2.22  0.9590  17/16.98  0.4400  02/  1.64  0.381D  16/15.58 

1.7783  3J.26  0.124D  25  0.7710  18/17.89  0.3050  02/  1.48  0  1870  17/16.27  0.8010  01/  0.90  0.6990  15/14.84 

J. 1623  37.22  0.1240  25  0.1650  18/17  22  0.548D  01/  0.74  0.353D  16/15.55  0.144D  01/  0.16  0.1260  15/14.10 

5.6234  41.19  0.1240  25  0.3320  17/16  52  0.981D  00/  0.01  0.648D  15/14.81  0.2580  00/  0.59  0.2270  14/13.36 

10.0000  45.14  0.1240  25  0.6410  16/15.81  0  1750  00/  0.76  0.117D  15/14.07  0.461D-01/  1.34  0.405D  13/12.61 

END  RINGBD  RUN 

STATEMENTS  EXECUTED'  1281233 
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Figure   2.      Prandtl-Meyer     Centered      Rarefaction        Fan        and        Breakdovm     Surface 

(Schematic) 
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STREAMLINE  2  R(M1)=54.4*A1 
BREAKDOWN  SURFACE  R(M1)=47.4*A1 
STREAMLINE  1         R(M1)=47.4*A1 

A1  =  (21/27rD2N0)"1[l*(y=l)M12/2]1/(y1) 


M^4 

y=1.54 


SUPPLEMENTARY   BREAKDOWN   SURFACE 


Figure   4.      Prandtl-Meyer        Flow       Field        Near      the     Corner, 
Streamlines    and   Breakdown   Surface 
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Including      Actual 
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FLUX  OF  SPECIES     H  [h]  =  .091 

M^4  y=1.54 

T0  =  2300(K)       WA  =  7.27  .        B  =  .05 
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DISTANCE  FROM  NOZZLE  LIP    Xs  (m) 
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Figure  5.   Flux  of  Species  H  at  Various  Stagnation  Densities 
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FLUX  OF  SPECIES     HF  [hf]  =  .091 

M-]=4  )f=1.54 

T0  =  2300  (K)      WA  =  7.27  B  =  .05 


DISTANCE  FROM  NOZZLE  LIP    X.  (m) 

L , L , \ J_J 
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10 


Figure  6.   Flux  of  Species  HF  at  Various  Stagnation  Densities 
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FLUX  OF  SPECIES     H2  [h2]=.104 

M-,=4  ^=1.54 

T0  =  2300  (K)      WA  =  7.27  B  =  .05 


DISTANCE  FROM  NOZZLE  LIP    X5  (m) 
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Figure  7.   Flux  of  Species  H2  at  Various  Stagnation  Densities 
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FLUX  OF  SPECIES     DF  [df]=.135 

M-^4  ^=1.54 

T0  =  2300(K)      WA  =  7.27  B=.05 


DISTANCE  FROM  NOZZLE  LIP    X§  (m) 
\  \  \         S 
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10 


Figure  8.   Flux  of  Species  DF  at  Various  Stagnation  Densities 
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FLUX  OF  SPECIES     He  [He]  =  .579 

M1=4  }C=1.54 

T0  =  2300(K)      WA  =  7.27  B  =  .05 


j£=.0075(kg  m-3) 


DISTANCE  FROM  NOZZLE  LIP    X<-  (m) 
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Figure  9.   Flux  of  Species  He  at  Various  Stagnation  Densities 
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FLUX  OF  SPECIES     HF  [hf]=.091 

M1=4  ^=1.54 

T0  =  2300  (K)      WA  =  7.27  B  =  .05 


DISTANCE  FROM  NOZZLE  LIP    X<-  (m 
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Figure  10.   Flux  of  Species  HF  at  Various  Values  of  Breakdown  Parameter 
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Figure  11.   Flux  of  All  Species  at  Typical  Operating  Conditions 
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Figure  12.   Schematic  Display  of   Complete   Breakdown   Surface   in  a  Ringjet 
Exhaust  Plume 
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